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Abstract. We compute the distribution of the partition functions for a class of one-dimensional 
Random Energy Models (REM) with logarithmically correlated random potential, above and at the 
glass transition temperature. The random potential sequences represent various versions of the 1/f 
noise generated by sampling the two-dimensional Gaussian Free Field (2dGFF) along various planar 
curves. Our method extends the recent analysis of [13] from the circular case to an interval and is 
based on an analytical continuation of the Selbcrg integral. In particular, we unveil a duality relation 
satisfied by the suitable generating function of free energy cumulants in the high-temperature phase. 
It reinforces the freezing scenario hypothesis for that generating function, from which we derive 
the distribution of extrema for the 2dGFF on the [0, 1] interval. We provide numerical checks of 
the circular and the interval case and discuss universality and various extensions. Relevance to the 
distribution of length of a segment in Liouvillc quantum gravity is noted. 
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1. Introduction 

Describing the detailed statistics of the extrema of M random variables Vi with logarithmic correlation 
built from those of the two-dimensional Gaussian Free Field (2dGFF) V(x), is a hard and still mostly 
open problem. It arises in many fields from physics and mathematics to finance. The 2dGFF is a 
fundamental object intimately related to conformal field theory [I], and being also a building block of 
the Liouville random measures e v ^ z,z ' dzdz attracted much interest in high-energy physics, quantum 
gravity, and pure mathematics communities, see [2] for an extensive list of references. In the context 
of condensed matter physics the 2dGFF is of interest to describe e.g. fluctuating interfaces between 
phases e.g. their confinement properties, multi-fractal properties of wave functions of Dirac particle 
in random magnetic field [4] and associated Boltzmann-Gibbs measures j5] , glass transitions of random 
energy models with logarithmic correlated energies [6], 2d self-gravitating systems [7] etc.. Descriptions 
of the level lines of the GFF as Schramm-Loewner Evolutions (SLE) and conjectured relations to 
the welding problem [5] have also contributed in revival of interest in the statistics of the GFF. In 
mathematical finance there is a strong present interest in limit lognormal multifractal processes [3j (also 
called log-infinitely divisible multifractal random measures), which is but a closely related incarnation 
of the same object, see e.g. [10j[TT]. Last, but not least important, is to look at the logarithmically 
correlated random sequences as those representing various instances of l/f noises, see e.g. [T2] and 
[13] • Such noises regularly appear in many applications, and were recently discussed in the context 
of quantum chaos, where logarithmic correlations arise in sequences of energy levels |14j or, as one 
can surmise, in the zeroes of the zeta Riemann function. All this makes understanding extreme value 
statistics of such noises an interesting and important problem. 

While the leading behavior V m i n ~ — 2AlnM is rigorously proved [15], surprisingly little 
knowledge exists on finer properties of the statistics of the GFF-related minima, even heuristically. 
To serve this as well as many other purposes it is of high interest to study the canonical partition 
function Z(fS) = X)f=i e~^ Vi for the corresponding Random Energy Model (REM) as a function of the 
inverse temperature f3 = 1/T. The distribution P(F) of the free energy F = —ThiZ reduces in the 
limit of zero temperature T = to the distribution of the minimum V m i n . A few instances of REM 
can be solved explicitly, and are frequently useful as approximations: (i) uncorrelated energies with 
variance ~ InM, i.e. Derrida's original REM [IB], which gives the correct constant A [l](ii) paths with 
random weights on trees, whose energies exhibit a similar logarithmic scaling of correlations, but with 
a hierarchical structure rather than a translationally invariant one [T71 [TS] (iii) the infinite-dimensional 
Euclidean version of logarithmically correlated REM and its further ramifications 0119]. In particular, 
the close analogy of GFF-related statistical mechanics with the models on trees [H [5J , also noted in 
probability theory [T5], arises naturally in an approximate, i.e. one loop, RG method, and led to the 
conjecture [5J that: 

V Tn i n = a M + b M y (1) 

with 

a M = A(-21nM + 7lnbiM + 0(l)) , b M = A+ 0(l/ln(M)) (2) 

where 7 = 3/2 and y is a random variable of order unity whose probability density has universal tails 
p{y) ~ \v\ eV on the side y — > — oo. In addition it was convincingly demonstrated that the log-correlated 
REM exhibits a freezing transition to a glass phase dominated by a few minima, at the same T c as 
predicted by (i) and (ii) [6] [19]. An outstanding problem left fully open was to characterize the shape of 
the distribution of the minimum beyond the tail, and in particular investigate whether the universality 
also extends to that regime. 

To address this issue, Fyodorov and Bouchaud |13j (FB) recently considered a particular circular- 
log variant of REM. Denoting here and henceforth the averaging over the random potential with the 
overbar, the circular-log model is defined via the correlation matrix Cy = ViVj identical to those of 
M equidistant points Zj = exp(z^-) on a circle Cjk = lG[Zj — Zfe), where G(z — z') = — In \z — z'\ is 
the full plane Green function of the 2dGFF. Equivalently, the above covariance function represents a 
27r-periodic real-valued Gaussian random process V(x) — YluLi \ v l e * fa + e_ ) with a self-similar 
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spectrum (viv m ) = l~^ 2H+1 ^ 5i m characterised by the particular choice of the Hurst exponent H = 0. 
Such a process therefore represents a version of the so-called 1/f noise. 

From the moments Z n FB reconstructed the distribution P(Z) above and at T c . From such a 
point they proceeded by assuming that for such a model the same freezing scenario as found in Ref. 
[6] holds so that the generating function 



gf3 (y)=eM-e^Z/Z e ), Z e = M 1+/3 /r(l — (3 2 ) (3) 

remains in the thermodynamic limit M 3> 1 temperature independent everywhere in the glass phase 
T < T c . As a result of such a conjecture they arrived at the distribution of the minimum of the random 
potential in their problem. The corresponding probability density for the variable y (defined in (|T|) 
with A — 1) turned out to be given by p(y) — —g'^iy) where 

goo (y) = gpM = ^ y/2 Ki(2ey/ 2 ). (4) 

Such a density does indeed exhibit the universal Carpentier-Le Doussal tail p(y — > — oo) ~ —ye y . 

Our broad aim is to investigate analytically and numerically the validity and universality of the 
above result, and to extend it to other models with logarithmic correlations. In pursuing this goal 
we will be able, in particular, to extract statistics of the extrema of the (full plane) GFF sampled 
along an interval, [0, 1], with eventually some charges at the endpoints of the interval. This breaks 
the circular symmetry of the correlation matrix and one finds a different distribution. The moments 
Z n turn out to be given in some range of positive integer n by celebrated Selberg integrals [20 [§] and 
a first (non-trivial) task is to analytically continue them to arbitrary n. After suggesting a certain 
method for such a continuation we are able to deduce the distribution of free energy P(F) and gp(y) 
at the freezing temperature f3 — [3 C . The same conjecture as in FB then yields the distribution of 
the minimum. As a by-product of our method we reveal a remarkable duality property enjoyed in the 
high-temperature phase by the generating function precisely defined as in ^ and unnoticed in 1 13] . 
We conjecture such a duality to be intimately related to the mechanisms behind freezing phenomenon. 
Finally we use direct numerical simulations to verify the freezing scenario for the circular ensemble 
and the resulting distribution Q, as well as to test the new results of this paper for the interval case. 
Universality and other cases are discussed at the end. 

2. Model and moments 

2.1. Interval model 

Our starting point is the following continuum version of the partition function of the Random Energy 
Model generated by a Gaussian-distributed logarithmically-correlated random potential V(x) defined 
on the interval [0, 1]: 

ri 

Z = e (i2 / dxx a {l-x) b e~ fiv{x) (5) 



o 



with a, b > —1 real numbers and j3 > 0. The potential V{x) is considered to have zero mean and 
covariance inherited from the two-dimensional GFF: 



V{x)V{x') = C(x - x 1 ) = -2 In \x - x'\ . (6) 

For the integral (|5| to be well defined one needs to define a short scale cutoff e <C 1. We therefore 
tacitly assume in the expression ^ V — > V e , with the regularized potential being also Gaussian with 
a covariance function C e (x — x'), such that the variance is C £ (0) = 21n(l/e). We put for convenience 
the factor e in front of the integral to ensures that the integer moments Z n are e— independent in the 
high-temperature phase, see Eq.([7]) below. At this stage we do not need to specify the e— regularized 
form [JJ] but it is convenient for our purpose below to require that C e (x) — C(x) for |i| > e. Note 

§ In a somewhat different but related context this fact was noticed, but not much exploited in [9] 

|| there are various useful cutoffs, e.g. the circle average, see e.g. [2], or the scale invariant cone construction, see e.g. 

m 
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that for a = b = the Gibbs measure of the disordered system identifies with the random Liouville 
measure, and that Z can be interpreted as the (fluctuating) length of a segment in Liouville quantum 
gravity see e.g. [2] . 

Below we will also consider a grid of M points Xi, uniformly spaced w.r.t the length clement 
dl = dxx a (l — x) b and the set of values Vi — V(xi), i = 1, ..M. The correlation matrix ViVj — Cij at 
these grid values are C v} = -21n(|i - j\/M) for i ^ j, and C u = 21nM + W where W = ln(l/(eM)) 
is a constant of order unity, and we will be interested in the limit 0j] of large M at fixed eM. 
This generalizes the grid on the unit circle studied in [T3] where Xj = e l0j with Oj = 2irj /M and 
Cij — C(xi — Xj) = — 21n(2| sin 2 |). We will compare below the two situations. In each case one 
defines the corresponding (discretized) REM by the partition function Zm — YlfLi e - ^. We expect, 
as shown in [TS] and discussed below, that there is a sense in which universal features of the discretized 
version are described by the continuum one in the large M limit. 



2.2. positive moments 

Let us now compute the positive integer moments of Z. Denoting 7 = (3 2 , a straightforward calculation 
gives 

pi pi n -. 

t—1 l<i<j<n J 

where the small scale cutoff is implicit and modifies the expressions for \xi — Xj\ < e. For a fixed 
n = 1, 2, .., a well defined and universal e — * limit exists whenever the integral Q is convergent, in 
which case it is given by the famous Selberg integral formula |20j Z n = s n , with: 

s " (7 ' a ' 6) = }J r[2 + a + 6-(n + J --2) 7 ]r(i- 7 ) (8) 

where T(x) is the Euler gamma- function. For a, b > the domain of convergence is given by 7 < \/n. It 
corresponds to the well known fact that for continuum REM models the distribution of P(Z) develops 
algebraic tails [+] hence integer moments Z n become infinite at a series of transition temperatures 
tJ"^ = s/n. The true transition in the full Gibbs measure happens however only at T c = 1 i.e. 
7 = 7 C — 1. Above T c the distribution P(Z) exists in the limit e = 0, while the formally divergent 
moments start depending on the cut-off parameter e. Analogous result arises in the log-circular 
ensemble [T3] where the moments of Zm were analyzed, as recalled below. The generalizations for 
complex a, 6, f3, which connect to sine-Gordon physics, as well as a detailed study of the competition 
with binding transitions to the edges for a, b < —1 (in presence of a cutoff) is mostly l eft fo r future 
studies, although some remarks about the binding transitions are made below in Section 3.2p" 



2.3. negative moments 

Our first aim is to reconstruct the distribution P(Z) from its moments in the high temperature phase 
7 < 1. This entails analytical continuation of the Selberg integral which is a well known difficult 
problem. Here we present a solution of this problem at T Cl the most interesting point. Let us first 
obtain the negative integer moments for any T > T c . It is convenient to define: 

z = r(l - i)Z = e~ 0f , z n = ^ (9) 

which, as found below, and in [T3], has a well defined limit as T — > T c + . One then checks for a = b = 
the following recursion relation: 

z n _ r[l - T17] r 2 [i - (n - 1)7] r[2 - (n - 2)7] 



z n -i r[2 - (2n - 3)7] T[2 - (2n - 2)7] 

% in practice we want that mini(|a;i — ajj+il) > e. 

+ for finite grid 1/M these tails are cut far away by log-normal behaviour, see a detailed discussion in [13] 

* the full conditions for convergence in (fSp are 5R(a), 5R(6) > — 1, < min(l/n, (a + l)/(n — 1), (b + !)/(" — 1)) 



(10) 
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with z\ — r(l — 7) (which also implies z = 1), and a similar formula for any a, b. Let us now perform 



the formal analytic continuation to negative integer moments mk = Z-k in the above recursion ( 10 1 
as mfe/rnfe+i = z„/z„_i|„^_fe. It is then easy to solve the recursion starting from too = Zq = 1. 
Restoring a, b we find: 

_ A r[2 + fl + b+(fc + , ? + i) 7 ] 

k / = i i r[i + (i-i) 7 ]r[i + a + i 7 ]r[i + 6 + j 7 ] 1 1 

We have checked that these expressions satisfy the convexity property zP~ m z™~ n > z^ n for any 
integers n < to < p of arbitrary sign, which is a necessary condition for positivity of a probability. For 



a = b = the formula (11 1 was announced very recently in |llj as a rigorous consequence of certain 



recursion relations for Selberg integrals. 

Note that the domain in a, b where ( 11 ) remains well defined extends to a > — 1 — 7 , b > — 1 — 7 , 
a region larger than the naive expectation a,b > —1. This is a signature of the competition between 
binding to the edge and the random potential as discussed below. 

2.4- From moments to distribution: the circular case and duality in the high-temperature phase 

Let us recall for comparison the corresponding analysis for the circle [13 . There, the corresponding 
Dyson Coulomb gas integrals give z n = T{1 — wy), and such simple formula admits the natural 
continuation to negative moments n — —k. This allows to immediately and uniquely identify the 
distribution of 1 j z and leads to the probability densities: 

P(z)^p- 2 z- 1 ^ 2 - 1 exp(~z~ 1 ^ 2 ) , P(/) = /3- 1 exp(///3-e ///3 ) (12) 

The latter formula implies that the free energy is distributed with a Gumbel probability density for all 
T > T c . Alternatively the (formal) series for positive moments gp(y) := e~ ze/3y = X^S^Lo ^"-V z n e nl3y 
is directly summed using T(z) = e~~ t t z ~ 1 dt into the following generating function 



/>oo 

gfi {y)= / dtexpi-t-ePyt-P 2 } (13) 
Jo 

What went unnoticed in |13j was the remarkable duality relation satisfied by the exact expression for 
this functiorjU 

9p(y) = 9i/p(y) (14) 
To see this directly define r = e^t -132 implying t — r t" 2 e~ y ^ , and after substituting this back to 



the integral ( 13 1 we see that 







/ dr 




>o 





1 f°° 1 1 

9p(y)=-M dTTl^e yl& vxv{-T-e ylf) T~lP} (15) 
P Jo 

exp{-T-e y/t3 T~^} = gi(y) (16) 

as second term in the integrand gives no contribution being full derivative of the expression vanishing 
at the boundaries of the integration region. This transformation is formal in the sense that the 
function gi/p(z) defined above for < 1 has nothing to do with the true generating function in the 
low temperature phase j3 > 1. Rather, it is just obtained by taking the formula valid in the high 
temperature phase and changing everywhere /3 —* 1//3. However the duality relation still gives a 
precious information, e.g. it implies that an infinite set of derivatives (fldp)™ gp(y) = for any n > 1 
odd at the self-dual point (3 — 1~. In particular the exact result: 

dp9p{.v)\ p=K =0 , for ally (17) 

shows that the "flow" of this function as a function of temperature vanishes at the critical point, quite 
consistent with a freezing of the whole function (with continuous temperature derivatives) . It is in fact 

() In general such duality holds for the transformation /3 — > 0^/(3 but we specialized in this paper to /3 C = 1 
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quite amazing that precisely this generating function gp(y) — exp(— e^ y z), with precisely this built-in 
temperature dependence, is both conjectured to freeze and shown to be self-dual. It is thus tempting 
to conjecture that freezing and duality are related, i.e. it is gp(y) and no other variation of it (such 
as e.g. replacing e@ v by any other function of both y and j3) which freezes because it is self-dual in 
the whole high temperature phase. The same type of self-duality relation, as we demonstrate below, 
extends to the interval case supporting the conjecture. 

Unfortunately, the direct methods of resummation which work for the circular case fail for the 
more complicated problem at hand, the interval [0, 1]. For this reason one needs to develop a more 
general procedure, which is done below. 



2.5. From moments to distribution: generalities 



Instead here we now define the generic moments Mp(s) = z 1_s , Mg(l) = 1 for any complex s, at 
fixed inverse temperature (3. In particular, the generating function of the cumulants for the free energy 
/ = — lnz is related to Mp(s) via 

oo n 

^^ n K = lnM^(l + S ) (18) 

n=0 

Definition of the probability density P(z) implies the relation 

» + 00 

e 2t P{e t )e- st dt= Mp(s) (19) 
which can be inverted as the contour integral: 

e- 2 *P(e-*) = — f e- st M (s)ds (20) 

2l7T J 

e.g. along a contour parallel to the imaginary axis s = So + iui, provided the integral is convergent, so 
being chosen larger than any singularity of the integrand. 

Further using the definition ^ the function gp(y) is found to satisfy the identities 

/+oo 
e^ s -^g {y)dy = M {s)T{s-l) (21) 
-oo 

g (y)=p-i e ey± J e-^M^T^-^ds (22) 

Hence once we know Afg(s) we can retrieve all the interesting distributions. Moreover, relation (21 1 
defines after integration by parts the generating function of the cumulants for the probability density 
defined by p (y) = -g' p {y): 

°° n /'CO 

J2 S ~W C = In / pp(y) e ys dy = lnM^l + |) + lnr(l + |) (23) 

1 n - J— oo P P 

Comparison with (JT8J) yields after recalling the series expansion for lnT(l + s) in terms of the Euler 
constant and Ricmann zeta-function C,{n) the following model-independent relations: 

V = J-7eT, F rC U>2=K+ (-!)"(«- l)!C(n)T", (24) 

This relation is valid at all temperature and comes only from the definition of gp(y). It is most useful 
at p = flc = 1, if we accept the freezing scenario. Given that in that case the l.h.s. freezes at its value 
at P = 1 then we easily retrieve all cumulants of the free energy for all T <T C just from the knowledge 
of 3/3=1 (y). Conversely, it is useful to test the freezing hypothesis in numerics, as we will see below. 

Let us now discuss how these moment relations reflect duality for the circular case. In the latter 
model Mp(s) = T(l + (s — 1)7), hence from (23 1 one finds: 

00 n 

V ^^ = 11^(1 + S p) + lnT(l+ 4) (25) 
' n! p 
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which is manifestly invariant by the formal transformation (3 — > 1//3. The latter fact implies, via (23 1, 
the self-duality for pp(y), hence for gp(y). Such an indirect method of proving self-duality for gp(y) 
has advantage when direct verification is difficult in view of cumbersome and/or implicit form for the 
generating function in the whole high-temperature phase. We shall see later on that it indeed works 
for the interval case. 

3. Analytical continuation at the critical temperature and distribution of minima on the 
interval 

3.1. no edge charges 

Let us keep focussing on the critical temperature p = 1. Denoting Mp = x[s) = M(s), gp = \(y) = g{y), 
we start with a = b — case (no charges at the end of the interval) for the sake of simplicity. 



For negative integer values s = 1 — n one finds from ( 10 1 after exploiting the doubling identity 
T(2z) = 2 2z - 1 r(z)r(l/2 + z)/*fw the relation: 

To continue this formula to any s we will use the Barnes function, which under some mild conditions 
is the only solution [5T] of: 

G(s + 1) = G(s)T(s) (27) 
with G(l) = 1. The Barnes function G(s) is meromorphic in the complex plane and has zeroes at all 
negative integers s = 0, —1, —2, ... [21]. It can be computed as: 

G(z) = (27r) (z_1) / 2 e _ 5( 2 - 1 )( z - 2 )+/o !_1 dxxip(x) ^8) 
where ip(x) = T'(x)/T(x), the integral being on any contour not crossing the real negative axis. Using 



(27 1 one finds the following analytical continuation for the moments, which is one of the main result 



of this paper: 



M(s) 



1 



Gjs + I) 
G(s) 



(29) 



G(5/2) 2 7r s - 1 r(s)r(s + 2) 

with G(5/2) = A _3 / 2 7r 3 / 4 e 1 / 8 2~ 23 / 24 where A is Glaisher-Kinkelin constant A = e i/i2-C'(-i) = 
1.28242712. To guarantee that this is the correct continuation, we have checked (i) positivity: M(s) 
given above is finite and positive on the interval s G [0, +oo[ i.e. all real mome nts n = 1 — s < 1 exist, 
(ii) convexity: on this interval d 2 In M(s) > (iii) convergence of the integrals ( 20|22 1 for sq > 1. The 



latter can be used to compute 3/3=1(2/) and P(f) = eJ /(2iri) J e~? s M{s)ds, which are plotted in FigJTJ 
Note finally that it reproduces the negative integer moment s (p~T|) for a= 6 = 0, 7=1. 

The free energy cumulants are to be determined from ( |23j24| , and one finds < y >= \ — 2je — 



ln(27r), < y 2 > c = ^| ^, and for general n > 3: 

< y n > c = (-)"- x (n - l)!(C(n - 1)(2" - 4) - CW(2"3 - k) + 2" +1 - 1 - 2" n ) (30) 
with k = 4 and the same formula for < /" > c with k = 3. As a comparison for the circle M(s) = T(s), 
hence < y >= 2 < / -2j E , and < y n > c = 2 < /" > c = 2(-l) n (n - l)!C(n) for n > 2. 

An important property of g(y) at criticality is its decay at y — > —00. Deforming the integration 



contour in (22 1 one obtains g(y) as a sum of residues over the (multiple) poles of M(s) at s = —n, 



which generates the expansion in powers of e v . 



g(y) = l+(y + A')e^> + (A + By + Cy 2 + ^y 3 )e 2 * (31) 
+ e vJ2 - i —d 2n e- s y2- ( - n+1 K 2n+3+is K n + 1 (32) 

n=2 ( 2n >- 

r(n + 1 + s) 2n+1 T(n + 3 + s)G(s + |) 2 M(n + 1 + s) 

X 



(s - l)s 2 (s + l) 4 ..(s + n - l) 2 "-!G(n + 1 + s + |) 2 
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Figure 1. Color online. Analytical predictions for the interval [0, 1] with no edge charge: (i) 
Left: plot of gp c (y) which, according to the freezing scenario is also, up to a shift, the cumulative 
distribution of the minimum V m i„ (ii) Right: the free energy density P(f) at the critical temperature 
C for the interval. Both are obtained by the appropriate inverse Laplace transforms \20\ and \1'2\ 
from the analytical continuation l|29| of the moments as indicated in the text 



with A' = 2j E + ln(27r) - 1 and C = -0.253846, B = 1.25388, A = -5.09728. Let us recall that for 
the circular model the expression Q implies: 

g^rc) {y) = l + e y {y _ l + 2 7b) + e 2t/ ( l y _ 5 + + (33) 

The behaviour g(y) — 1 ~ ye y is precisely the universal tail found by Carpentier and Le Doussal [5J. 
It has its origin in the 1/z 2 forward tail which the probability density of z develops at critical (3=1, 
with the first moment < z > becoming infinite. On the other side y — > +oo one expects much faster 
decay, for example g^ arc ^(y) = v^e* ~ 2e " /2 (1 + ^e _y / 2 + ..). 

3.2. extension to edge charges, binding transition 

Extending these considerations for any a, 6, one finds: 

Mis) = 22s 2 +s(i+2(a+W-3-2(a +j , V - s G(2 + a)G(2 + b)G(4 + a + b) 
U r(2 + 2±*)G(2 + s±k) 2 G(l + 

r(i + + s )G(i + z$p + S ) 2 G(| + + s ) 2 

X G(s)G(l + a + s)G(l + b + s)G(3 + a + b + s) ( ' 

and checks again positivity and convexity for s G [0, +oo[ (for a,b > —1). We give only: 

— c 7T 2 

y 2 a ,b = y + IE + 3^(4 + a + b) - 0(2 + a) - 0(2 + b) (35) 
with </>(a;) = ')p(x) + (x — l)^'(x), and the case a — b in the limit a — > +cxd where one then finds: 

^; a =ln(8a) + ^+ 7 +l + 0(a- 1 ) 

K,a = -^ "2C(3) + 0(a- 1 ) + .. (36) 

i.e. all cumulants have a limit except the second one. This limit is discussed again below. 
A remarkable case is a = 6 = —1/2. Then a simplification occurs: 

M(s) = M_ 1/2 ,_ 1/2 ( S ) = 2 2s2 - s -V"*^i^ (37) 

One can trace this simplification to the fact that the structure of the correlation matrix becomes much 
simpler in that detailed in Appendix A. The corresponding distribution is easlily found as 
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(again this is for (3 = j3 c = 1): 



/ \ /71\c./9 1 [ Z dzi [ + °° dt _*2__q«/hT2*_«; J_ p -2vTIT2t 

p W-<8» 3 ' 2 F(wi„ w*' ' (38) 



"1 

.,2 



which reproduces the above moments, and behaves as P(z) ~ 7r/z at large z. This yields, after some 
manipulations: 

TT f +0 ° dt o./uTo, Z" 00 A e» 



ff (y) = - / « e-^-a^t / ^ - fl e -Vw5«-^ dw (39 ) 

and one finds = 41n2 - 3 + 2tt 2 /3 = 6.35232 and y = 1 - 2 7i; - ln(f ) = -0.606014 hence 
/2 C = 41n2-3 + 7r 2 /2. 

Let us now discuss briefly the case a, b < —1, for simplicity we focus on b — a. In that case, 
the model ^ requires, at least naively, a short scale cutoff to avoid the divergence near the edges. 
However, from e.g. the discussion of Appendix D in [B] we know that there should be a competition 
between the random potential in the bulk and the binding effect by the edge: in presence of disorder 
it may be more favorable for the particle to explore the bulk and to remain unbound from the edge. It 
is quite nice that our analytical continuation captures that effect. As mentioned in Section |2.3[ from 
the negative moments one can guess that the complete domain over which the high temperature phase 
extends is: 

a > -1 - 7 and 7 < 1 (40) 

with 7 = /3 2 (here (3 C — 1), where equality in the first condition corresponds to the binding transition 
to the edge, while the second to the freezing transition. This implies in particular that for 7=1, the 
case studied above, the binding transition occurs at a = —2, and that for any larger value of a the 
system should be at bulk critical freezing, with however some continuous dependence in a. We can 



indeed check that the result (34 1 for M(s) leads to a well defined probability P(z) for any a > —2. 



For instance one sees that the formula (35) yields a finite y 2 for any a > —2, which however diverges 
as a — > — 2 + . The domain of definition becomes s > — 1 — a for — 1 > a > —2, as the resulting P(z) 
acquires now a broader tail ~ l/z 3+a at large z, while it was ~ 1/z 2 for a > — 1. As a — > —2 the 
tail becomes non-normalizable as ~ 1/z, a signature of the binding transition. The case a = —3/2 



provides a good illustration as ( 34 1 again simplifies into 



M(s) = M_ 3/2l _3/ 2 (s) - 2 2s2 - 5s+ V/ 2 - s r( S - 1) (41) 

which implies that the random variable z can be written z — zie~^ 2 where Z\ > and fi are two 
independent random variables, fa being gaussian distributed with f% = — ln(27r) and /f = 4 In 2, and 
Z\ with distribution P\{z\) — z 1 3 ^ 2 e _1//zi /y / 7r, leading to the explicit form: 

1 Z" 00 ,31, (t + ln(27r)) 2 N , . 

V ; z^TrVSh^J-^ 2 z 81n2 ; v ' 

which does exhibit the ~ 1/z 3 / 2 tail at large z. We leave further studies of the global phase diagram 
for arbitrary a, b to the future. 

4. High temperature phase for [0, 1] interval with no end charges 

Let us consider the segment [0,1] at any j3 < f3 c = 1, i.e. 7 = f3 2 < 1. The moments must satisfy 
(using again the doubling identity): 

m ( s + 1) = 2 2 +^ r(§ + S7 )r(i + I + S7 )r(f + \ + SJ ) 

Mp(s) tt r(l- 7 + s 7 )r(2 + 7 + s7)r(l + s 7 ) 1 ' 
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We need to find a way of continuing the moments to the complex plane. To this end we define the 
function Gp{x) for SR(a;) > by [23]: 

Q/2^ /0 _, , f°°dt, e-^-e'^ , e~* 2 Q/2 



ln W ^ln(2,) + j o T( (l-e^)(l-e-^) + ^ 2 -^ + ^) 
where Q = + 1//9- This function is self-dual: 

= (45) 
and satisfies the property that we need, see e.g. [23] and Appendix B, 

G (x + 0) = p 1 l 2 -^(2^T(Px) G {x) (46) 

One can check that Gp{x) for j3 — (3 C — 1 coincides with the Barnes function G(x) defined in the 
previous Section, e.g. setting (3 = 1 in (46 1 one sees that Gi(x + 1) = T(x)Gi(x), and, using Q = 2 we 
have Gi(l) = 1. Similarly to the standard Barnes function the new function Gp(x) has no poles and 
only zeroes, and these are located at x — ~n[3 — m/ (3, n, m = 0, 1, ... It provides us with a natural 
generalization which can be used to perform the required analytical continuation for any temperature. 
Using the above properties we find that 

M M- A gC-DC^^DV- Eg + 1))G ^ + I + PS)G ^ + ^ ( 2 tklM 

M) " G^+l + MGj^ + Ps) 

(47) 



with 



G a + ffG {2(3+l) 



Ap G?(f + I)G^ +/ 3) G/3 (f + ^)) 

reproduces correctly the recursion relation ( 43 1 , hence provides an analytical continuation for the 
moments valid for j3 < f3 c = 1. We have checked numerically that it does satisfy positivity, convexity 
and a convergent inverse Laplace transform from which one can compute P(z) and gp{y) using (20 1 
(22 1. We will not study these in details here, but give only a few properties. 
Let us first check the duality. One easily sees that if one defines 

Mp{s) = 2 1 ~ s M (s) (49) 

then hxMp{l + §) + hiL(l + ■§) is fully invariant under (3 -> 1/(3. From (23 1 it implies that all y" c 
with n > 2 are invariant by duality, only the average y is not. This is not a problem since this average 
is not expected to be universal, and is easily remedied by defining z = z/2 (which could have been 
done from the start) and gp(y) = exp(—e l3y z/2). Hence we conclude that up to such a trivial shift 
the probability p,a(y) = —g'p{y) is self dual, i.e. ft\/p{y) — pp(y). From the discussion in the previous 
Section we conjecture that it is this function which freezes at (3 = (3 C = 1. 



From the result (|47j) we can extract the cumulants of the free energy using ( 18 1 . We only discuss 
here the lowest non-trivial cumulant, given by 

T = ¥ C - yT 2 = ^A 2 ln M p (l + s)\ s=0 (50) 

3/? 1 3 3/3 3 1 9 R 2 tt 2 

= 41n2 + h (^ + -) + h ([3+-) + h p (f + -)-2h (P+-)-h p (2p+-) + ^- 

(51) 

where we have defined the self-dual function (see Appendix B): 



DC- 



dt . xt A t 2 



hp(x) = h 1/p (x) = %\n.Gp{x) =Ux + J -e~ xt \l - _ _ & \ (52) 

and we have used ^'(1) = t 2 /6. The resulting curve f 2 as a function of (3 is plotted in Fig. [IT] One 
finds that it increases from J 2 "^ -> 0) = 3 to /2 c (/3 = 1) = 7tt 2 /6 - 27/4 = 4.76454. More discussion 
is given in Section and Appendix C, together with high temperature expansions. 
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5. Gaussian weight model 

We now briefly discuss a case where the above considerations fail, and present below some hints of 
why this may happen. 

We consider now the continuum partition function for the log-correlated field on the full real axis 
but with a gaussian weight: 

1 



^ * / dx e~ x V2 e -/3V» 
v2tt J-oo 



This problem is appealing as it leads to Mehta integrals and moments z^ — z n — Z n T{\ — /3 2 )" = 
r[l — j(3 2 ], i-e. simpler expressions than for the interval case considered above. 
At criticality (3 = 1 this implies M^ G \s + \)/M^ G \s) = l/r(,s) for s = —n, which naturally 
suggests M( G \s) = 1/G(s). This is positive for s > but, surprisingly convexity fails for for 
s > s c = 1.92586... Hence this is not an acceptable analytic continuation. 

To get another handle on the problem one notes that this model can be obtained from the large 
a limit of the interval problem [0, l] aa - Writing x — 1/2 + y and performing the change of variable in 
Q one finds: 

lim (27r)-"/ 2 2 2a "(8 a )t- M ^^„(a,a, 7 ) = z I f)( 7 ) (53) 

a — >+oo 

Not surprisingly one finds that the pointwise limit: 

M i - G \s)= lim (2 7 r)-( 1 -^' /2 2 2Q ( 1 - s )(8a)5 (1 - s2) M a , a (s) (54) 

a — >-+oo 

yields 1/G(s) as expected. From this we also get that for large a: 

d 2 s In M^ G \s) = -ln(8a) + <9 2 In M a . a (s) (55) 

While the second term is nicely positive for all s > the additional factor — ln(8a) makes the total sum 
negative for s > s c , violating convexity. In other words while M ayCL (s) corresponds to a well defined 
distribution of probability, corresponding to the problem on the interval with edge charges, Mq{s) 
corresponds then to this probability " convoluted by a gaussian of negative variance" and fails to be a 
probability. Note that such a shift in the second cumulant y 2 is indeed needed to obtain a finite final 



result in (36 1. All higher cumulants y n \ aa with n > 3 have a nice finite limit as a — > oo, and can be 



extracted from the generating function 

J2 -iy 1+n = - ~ ( s - l )^ s ) + s - o ln ( 2?r ) - o 

nl s 2 2 

n=0 

obtained from l/G(s). Hence the main problem seems to lie in the second cumulant, and one may 
speculate that it is related to an inadequate treatment of zero mode fluctuations. Another (possibly 
related) observation is that for a >> 1 the whole contribution to the [0, l] aa integral comes from a very 
small vicinity (of the widths of L a ~ W a ) °f the mid-point x = 1/2 of the integration domain. One 
expects a competition between L a and the regularization scale for the logarithm, so it may be that 
the result depends on the order of limits e — > and a — > oo. We leave further study of this problem to 
the future and now turn to numerical studies. 

6. Numerical study 

6.1. circular ensemble 

We now turn to the numerical checks for the random variables Vi on i = 1, ..M grid points and their 
associated REM of partition function Zm — X)i=i e~@ Vi . We start with the log-circular ensemble and 
study the M X M cyclic correlation matrix (choosing here W = 0): 

Q J = -21n(2|sin ^ ( ' M j) |) i^j, C u = 2\nM+W (56) 
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whose eigenvalues Aj, = 21nM — 2 ^^J^ 1 cos{|jnfc} ln{2 sin -^n} are all positive, with the uniform 
mode Aq = for any M. Let us recall that the relation to the continuum model defined above 
was established in [T3] where it was shown that at large M one has Zfy = z n Z™ for f3 2 n < 1 and 
Z^j ~ M 1+n & for p 2 n > 1 (the positive moments which formally diverge in the continuum). 
The random variables Vi are generated (for M even) as 




where the Xk and yk are two uncorrelated sets of i.i.d. real unit centered Gaussian variables. This is 
done using Fast Fourier Transform (FFT). 




4 4.25 4.5 4.75 5 0.05 0.1 0.15 0.2 0.25 

log(4 7i log(M)) 1/Log(M) 



Figure 2. Color on li ne. Left: Finite size sc alin g of ajvf for variables with logarithmic correlations, 
circular ensemble Eq. \56\ , and interval Eq. \60\ , from M = 2 8 to M = 2 19 . The predicted slope is 
7 = 3/2, numerically we find 7 = 1.4 ± 0.1. This is compared with independent random variables, 
the standard uncorrelated REM, where the prediction is 7 = 1/2 as observed. Right: Finite size 
effect for fejvf for variables with the same correlations. The data are consistent with a convergence as 
1/ logM and extrapolate to b^j = 1 ± 0.02, consistent with the predicted value b^j = 1 in each case 
(which means an unrescaled variance V^ in in agreement with the prediction given in the text, in 
each case). 

From the distribution of the minimum V m in in systems of up to M — 2 19 we have computed the 
coefficients ajvf, bu and the distribution of the variable y in by fixing y and the variance y 2 to 
their value for the distribution Q. The asymptotics of the coefficients au and bu in Q are shown 
in Fig. [2] They exhibit a reasonable agreement with the conjecture Q with A = 1 but one clearly 
sees that convergence is slow. Convergence to &m = 1 would mean that the prediction Vji in = 7r 2 /3 
is correct. The cumulative distribution Qm(u) of the rescaled minimum, i.e. the variable y 7 is shown 
in Fig. [3] where the cumulative distribution Q has been substracted. One sees that although the 
difference is small its convergence, if any, to zero is extremely slow (empirically a ~ 1/VlnM seems to 
roughly account for the data, but we do not wish to make any strong claim here). 

Then we computed the distribution of the free energy at various temperatures. In Fig. [4] we have 
first normalized the free energy distribution to the same average and variance as the unit cumulative 
Gumbcl distribution, i.e. exp(— e x ), then plotted the difference between the resulting cumulative 
distribution Q re sc{f) and the Gumbel expression. This shows that the convergence is very fast at 
j3 = 1/2 but rather slow already at (3 — 1, where we have little doubt for the result. This is consistent 
with the fact that the convergence for the minimum is so slow. 

To test the freezing scenario we also compute numerically gp{y) for various temperatures. First in 
Fig. p^we test the convergence of the numerically determined gp c= \{y : M) to the analytical prediction 
gp e (y) in Q as a function of M. Then in Fig. [6 we test whether gp(y,M) — gp c (y,M) at fixed 
[3 > p c decreases to zero as M becomes large, which is the freezing conjecture. In practice we 
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Figure 3. Color online. Circular case: cumulative distribution of the rescaled minimum Qm(v) 
minus the prediction Q based on the freezing scenario gp a {y)- The number of samples is 10 7 . The 
difference is small compared on the scale of unity. Although it is slow, the convergence is apparent. 
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first compute the free energies /j, compute their mean / and variance a, define rescaled energies 



/■=(/-/ + IeT - 2 lE )< s J^{\ - T 2 /2)/ y /o z and define gp{y,M) as the mean of e~ e * which, 

by construction and virtue of (24 1, has then the same average, —27^, and variance, 7r 2 /3, as gp c (y) in 
Q. Comparing Fig. [3] and Fig. [5jwe see that a good fraction of the difference in Fig. [3] is already due 
to finite size corrections at (3 C (which have nothing to do with the testing the freezing scenario). 





Figure 6. Color online. Circular case: direct test of the freezing scenario: convergence of 
9/3(3/1 Af) — 98 c {Vi Af) (both are numerically measured and rescaled as explained in the text). We see 
that the scale is smaller by a factor around 2 to 4 than on Fig. [3] but that convergence is very slow. 
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6.2. universality of circular ensemble: cyclic matrices, GFF inside a disk with Dirichlet boundary 
condition 

It is important to discuss now the universality of this result, as it is a rather subtle point. The general 
issue of universality for logarithmic REM's can be formulated as follows. Consider sequences of M- 
dependent correlation matrices Cy\ What are the possible universality classes for the associated 
REM in the limit M — ► +00, what are their basin of attraction and conditions for convergence? 
One may ask two questions: (i) extremal universality classes, i.e. correlation matrices which 
have asymptotically the same distribution of the minimum V m i n (up to a shift by a M-dependent 
constant) (ii) more restrictive universality classes valid for any /3, i.e. correlation matrices which have 
asymptotically the same distribution of free energy, and generating function gp{y) (up to a shift by a 
Af -dependent constant) for any (3. It is reasonable to expect each latter universality class (ii) should 
correspond to a continuum model. Obviously, two sequences which belong to the same class 

(ii) also have the same distribution of extrema. But there are counter examples to the reverse (see 
below). Classifying these classes being a formidable problem, here we only make a few remarks about 
the universality class of the circular ensemble. The class corresponding to the interval is discussed 
below. 



Let us start from ( 56 1 and discuss various generalizations in the subset of cyclic (called also 



periodic or circulant) matrices, i.e which can be written as: 

Af-i 

fe=0 

with (Af-dependent) real eigenvalues Afc. The eigenvalue Ao corresponds to the uniform mode (often 
called zero mode in the GFF context). Logarithmic correlations mean that we assume that A^ ~ l/k 
in some broad range of k at large M, as specified below. 



Starting from (56 1 let us first make the observation that adding a fixed W > of O(l) on the 
diagonal of Cy shifts all eigenvalues by a constant 0(1) and does not change the universality class, in 
both sense (i) and (ii), at large M. On the other hand, a shift: 

Cy -> Cy + a (58) 

for all shifts only the the uniform mode Ao — > \q + o. It is equivalent to add a global random 
gaussian shift v to all Vi, i.e. Vi — > Vi + v where a = v 2 . It thus results in the the convolution of the 
distribution of V m in (and of the free energy) by a gaussian of variance a. One such example, discussed 
again below, is to consider the distribution of the GFF (using the full plane Green function) on a circle 
of radius R < 1 (and cutoff Re, i.e performing a global contraction): it shifts all Cy — > Cy — 21ni? 



in (56 1. Hence we keep in mind that there is really a family of distributions differing by their second 
cumulant, and will enforce in our numerics the condition Aq = which we believe selects the distribution 



6.2.1. GFF along an arbitrary circle One possible generalization of the circular model ( 56 1 along these 
lines is the GFF inside a disk of radius L with V — on the boundary as studied by e.g. Duplanticr 
and Sheffield [2J. Using the Dirichlet Green function Gl(z,z') = — In i^Z^Ji j the correlation matrix 
for the discrete model on a circle of radius R inside the disk is then for j and denoting p = R/L: 



2o\ sinf — j )\ 

C . j = -2hx— p \ 2 /' — , C, J = 21nL + 21n(l-p 2 )-21ne (59) 
V 1 + p — 2p cos(Oi — 0j) 



In the small p = R/L limit, equivalently fixed R and large L one finds Cy ~ — 2 In p — 21n(2| sin( ' ~ 3 ) \ 
and Cu « 21nL — 21ne. Choosing [26 e = R/M one sees that one recovers indeed the FB model (56 1 
(with W = 0) up to a shift a = 21n(L/i?) in the zero mode Ao of the matrix, i.e all eigenvalues of the 
correlation matrix are the same as FB except the uniform mode. This gives us the precise meaning 
of the universality of the results of FB [13 : it holds for small p = R/L for the Dirichlet GFF on the 
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Figure 7. Color online. Eigenvalues of the correlation matrix corresponding to the GFF along a 
circle in a disk domain with zero boundary condition (Dirichlet). 



disk and up to a (trivial) convolution by a gaussian of width 21n(L/i?). The next question is whether 
the universality extends to other circular contours on the disk with R/L not necessarily small. The 



answer is no, as can be argued from examination of the eigenvalues, diagonalizing (59) for arbitrary p. 
As shown in Fig. [7] at large M the eigenvalues are essentially the same as the ones for FB, i.e. small 
p, apart from the few largest ones - whose number does not change and remains finite as M becomes 
large. We expect that however that since these are the largest eigenvalues, despite being few they will 
change the distribution of the maximum which hence will depend continuously on the ratio R/L (with 
a similar discussion as above concerning the zero mode and convolution by a gaussian) . 



100 




Figure 8. Color online. Eigenvalues of the correlation matrix corresponding to the periodic models 
defined in the text. 



6.2.2. other periodic models On the other hand a much stronger universality property appears to 
hold when only the smallest eigenvalues are changed. Hence we now test whether the results obtained 
for the circular case remain valid for all periodic cases (57 1 with the same behaviour of A& ~ X/k. 
Again it is important that Ao be fixed to zero. If Ao > this amounts to convoluting the distribution 
of the minimum by a Gaussian of variance Ao- For the model to be logarithmic and strong universality 
to hold we require that Afc — * X/k as M — * oo for < k <C M. We have tested this conjecture for two 
models. 

Model 1: Sharp model (SM) A fe = M/k for k = l,...,M/2 and A fe = M/(M - k) for 
k = M/2, . . . , M - 1 



Statistical Mechanics of Logarithmic REM: Duality, Freezing and Extreme Value Statistics ofl/f Noises generated by Gaust 







CO. 



o 



008 
006 
004 
002 


002 
004 
006 



008 



7 



Circular M=2 
Model 1 

Model 2 \ 
Circular M=2 10 
Model 1 
Model 2 




-10 



Figure 9. Color online. Universality in the case of periodic (circulant) correlation matrices: 
cumulative distribution of the minimum, with subtraction as in Fig. [3] convergence to a common 
curve is faster than to the global analytic prediction. 



Model 2: The long range model (LRM) which is some discretization of the Joanny-deGennes 
elasticity of the contact line [53], = 2tt/ \/2(l — cos(27rfc/M)). 

The eigenvalues of these models are compared to the one of the circular case in Fig. [8] and one 
can see that they differ only for k near M/2. As can be seen in Fig. [9]the convergence of these models 
to the circular case at fixed M is much faster than their (common) convergence to the analytical 
prediction. We take this as a signature of the strong universality with respect to variations of the 
correlation matrix which change only the smallest eigenvalues, within the cyclic class. We check in 
the Appendix C that the first terms in the expansion of f 2 are the same for all these models which 
supports that the universality holds at any /3, i.e. both in the sense (i) and (ii) defined above. 



6. 3. interval 

We now discuss the [0, 1] ensemble. We take for the correlation matrix the Toeplitz form CV,- = C(i—j), 
i,j = l,..M: 

C7 4l = 4^(-l) fc log- + ^ <7 i?y = -21ogL--^ (60) 
i 

with W — 0. This matrix is not diagonal in Fourier space and we cannot use the FFT method. In 
practice we find the eigenvalues Afc and the normalized eigenvectors ipk{i) by a direct diagonalization of 
the matrix Cy. We then generate the correlated random potential as Vi = X^lc/ V^k~ x kfpk('i), where 
the Xk arc i.i.d. real unit centered Gaussian variables. Performing this sum together with the direct 
diagonalization is numerically expensive and limits the size of the number M of correlated numbers. 
In order to achieve a good statistics (~ 10 7 samples) we analyze here data only up to M = 2 12 . 



To justify our choice for the diagonal clement in ( 60 1 let us recall a useful property of any Toeplitz 
matrix: if the function f(9) — Cxi + 2 J2k=i Ci,fe+i cos (^0) ^ s positive V9 £ [0, 2tt), then C is 
positive definite (for any M). This is seen by noting that for any vector Vk, k = 0, ..M — 1, one 
has Elo VkVeC(k -l) = ^ ^f(9)\v(9)\ 2 where v{9) = E^" 1 v k e* ek and C(k) = ^ fj{9)e™ 



More importantly it can be shown that the reverse is true for large M [25]. For the choice of Eq(60l 



this function has a global minimum at 9 = 7r, for which f{6 = 7r) = 0. As a result the matrix CVj is 
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positive definite and in the large M limit, one finds that the smallest eigenvalue goes ra pidl y to zero 
and the eigenvector components alternate as (— l) fe . Note that the diagonal element in (60) behaves 
as Cu ~ 21nM + O(l) at large M, hence as expected, and similarly to the circular case. Though a 
convenient choice to prove positivity, there are other choices with similar behaviours at large M which 
would do as well. 

We have analyzed the distribution of the minimum V m i n and computed the coefficients clm, i>M an( l 



4 2_ 27 
4 



the distribution of the variable y in ill by fixing y = 7/2 — 2je~ ln(27r) and the variance y 2 = | 
to their value given by the analytically prediction. The convergence to 6m = 1, shown in Fig. (right), 
is thus a test of our prediction V^ in = |-7r 2 — 2 j. The convergence of the coefficients om and bu is 
quite similar to the circular case. The cumulative distribution Qu{y) of the rescaled minimum, i.e. 
the variable y, is shown in Fig 10 where the cumulative distribution of Fig[l](our analytical prediction) 
has been substracted. Again, the behaviour resembles the one for the circular case. 

The discussion of the universality for the interval class is more delicate since now the lowest 
eigenvector is no more generically the uniform mode. However a way to realize it from the GFF can be 
suggested similarly to the above discussion. One can consider the interval embedded near the center 
in a large disk with Dirichlet b.c. In the limit of small ratio p of interval size to disk radius the above 
interval model applies, again up to a convolution by a gaussian of variance 21n(l//?). 



CO. 



O 




Figure 10. Color online. Interval case: cumulative distribution of the rescaled minimum Qm(v) 
minus our analytical prediction, g/3 a (y), shown in Fig. Ill and based on the freezing scenario. The 
number of samples is 10 7 . The difference is small compared on the scale of unity. Although it is slow, 
the convergence is apparent. 



6.4- temperature dependence of the second cumulant of the free energy 

Finally we have also performed some numerical tests of the temperature dependence of our analytical 
results in the high temperature phase. We have computed numerically, and plotted in Fig. |ll| and [T2] 
as functions of (3, the variance of the free energy distribution f 2 as well as y 2 for the circular case 
(56 1 and f 2 for the interval case (60 1. They are compared to the analytical predictions, i.e. (i) for 
the circular case: 

F = (^ 2 /6)/3 2 (/3<1) and J 2 " = (tt 2 /6)(2 - T 2 ) (/?>!) ( 6 1) 
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which via Eq. (24 1 corresponds to y 2 — \{P 2 + 4s) for (3 < 1 which freezes into y 2 = n 2 /3 for 

(3 > 1, and, (ii) for the interval case formula H for f3 < 1 and ~p C {(3) = J 2 ' \/3 c = 1) + ^(1 - T 2 ) 
for (3 > 1. One can verify the good convergence in the high temperature phase. Questions related to 
the behavior for small (3, and how the numerical convergence could be further improved, is discussed 
in the Appendix C. 




Figure 12. Color online. Circular case: second cumulant y 2 as a function of inverse temperature 
f3 for various sizes M as compared to the analytic prediction given in the text. 



6.5. more open questions on universality 

Let us now indicating a simple example where universality (i) of distribution of minimum and of (ii) the 
free energy at any temperature, discussed above, may differ from each other. Consider the continuum 
problem on the circle but with an arbitrary smooth and a non singular weight < p\ < p(9) < p2 < 1: 

Z = e? 2 d0p(9)e-^ v ^ = ^ [** dBe'^^'^ (62) 
Jo Jo 

and we consider for instance 9 = f(9) = 6 + asin(0) with a < 1 and p(9) — 1/ f'(9). From the second 

form in (jSJ) one sees that the associated REM can be chosen as Zm = Yli z~^ Vi with correlation matrix 

Cij = — 21n |2 sm(^f(9i) — \ f{9j))\ for i ^ j and 9i = 2iri/M, neither a circulant nor Toeplitz matrix. 

As shown in the Appendix C, at small a, f 2 = \a 2 + 0(a 4 ) + 0((3 2 ), hence the free energy distribution 

clearly depends on a. On the other hand, the first form in (pi) suggests that it should have the same 
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distribution of the minimum V m i n as p(9) = 1. Indeed, V(8) — T\np(9) as T — > should have same 
extremal statistics for a given e regularization (with T — > before e — > 0) as V(0) provided p(6*) is non 
singular. The question of how the freezing scenario works under such a circumstances and what are 
the universality classes is left for future studies. 

Finally, another challenging question about universality, related to the GFF, is about REM's 
constructed along different curves in the plane than the circle, or the interval, i.e. V$ = V(z{) where 
the Zi lie along a curve and sample it at large M with a density described by some given arc length 
•\/dzdzp(z, z). One can use conformal maps to relate various curves to each others, e.g. a circle to a 
slightly deformed circle, with different weight functions p{z,z). Hence we are back to understanding 
the type of problem described in the preceding paragraph, and one should expect some universality in 
the distribution of the minimum. 

7. Conclusion 

To summarize, we have studied analytically and numerically random energy models based on Gaussian 
random potentials with logarithmic correlations. We have extended the Fyodorov-Bouchaud (FB) 
results from the circular ensemble to the interval. We have found the proper analytic continuation 
from the positive integer moments of the partition function, expressed as Selberg integrals, to arbitrary 
moments. This analytic continuation of the Selberg integrals, previously an outstanding open problem, 
is solved here. The solution involves Barnes functions and their generalizations which appear in studies 
of the Liouville field theory, hence strengthening the already noted link between the two problems. 
This solution, valid in the high temperature phase, allowed us to obtain the full distribution of the 
free energy / for (3 < f3 c and up to the critical point. It was generalized, at (3 — f3 c to the case where 
additional charges exist at the end of the interval. 

The knowledge of the generating function gp(y) — exp(— e 13 ^^^) at (3 — (3 C allowed us, via 
the same freezing scenario hypothesis as put forward in FB for the circular ensemble, to obtain the 
distribution of the minimum of the gaussian free field (GFF) on an interval, expressed as an integral 
transform of a Barnes function. The freezing scenario, which asserts that precisely this generating 
function gp(y) becomes temperature independent in the glass phase for (3 > (3 C , was until now based on 
a traveling wave analysis. While rigorous for the Cayley tree based REM, for which it was introduced, 
it was only based on a one-loop RG analysis for the type of models at hand [BJ. Here we made what 
we believe should be considered as a step towards better understanding of this freezing scenario: we 
discovered that, both for the circular ensemble and its interval counterpart, the analytic expression of 
gp(y) obeys in the high temperature phase the duality with respect to the transformation {3 — > 1/(3. It 
implies in particular dpg^{y) = 0, for all y, at /3 = [3 C in perfect agreement with a continuous freezing 
scenario. While one may notice that the generating function gp(y) is special as being the partition 
function of the Liouville model (see e.g. the discussion in [BJ) further connections to duality in Liouville 
field theory remains to be understood (the high temperature phase being the analogous of the weak 
coupling phase in Liouville). 

Detailed numerical calculations of the free energy distribution and of the function gp(y) associated 
to discrete REM versions of the circular and interval models were performed. The freezing scenario 
is consistent with our results, in both cases, though convergence is found to be very slow. The 
numerically obtained distribution of the minimum V m i n of M random Gaussian variables Vi with 
logarithmic correlation matrices CV, is found to lie close to the predictions, but with only very slow 
convergence as a function of M. In the high temperature phase the convergence to the FB result for 
the circular case and to the present one for the interval is found to be very convincing, and in full 
agreement with various high temperature expansions also performed here. 

The important question of the universality classes for discrete REM based on logarithmic matrices 
Cij and for their continuum analogs, is discussed. The continuum circular ensemble of FB is found to 
provide a single universality class for all circulant matrices with appropriate behavior of their spectrum 
at large M, for which we provide several examples. This strong version of universality holds for any 
temperature, i.e. identical distribution of free energy for all (3, up to a shift. A weaker version of 
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universality, holding only for the distribution of the minimum V rn i n is discussed through an example. 
As far as the connection to the GFF is concerned, we discuss the case where the field is sampled along 
a circle of radius R inside a disk of radius L with Dirichlct boundary condition. We demonstrate 
universality, up to the convolution by a gaussian, in the limit of a small ratio R/L, while in general 
the distributions depend on the aspect ratio R/L. 

The present progress opens many more fascinating questions. First one would want to extend these 
results to other curves in the plane, and even to two dimensional regions. The simplest extension, i.e. 
the case of the real axis with gaussian weight, also studied here, and for which the present methods 
are found to fail, shows that more remains to be understood before this can be achieved. Unbounded 
regions seem to pose a problem, and so does the control of the zero mode. The question of classifying 
the universality classes remains as a tantalizing open question. One can expect that the conformal 
invariance of the 2d GFF will play a crucial role in that classification, as it allows to map one curve 
into another one, with a change in the local length element. The question of which models obey duality 
and what is the precise connection to the freezing scenario is also outstanding. Further exploration of 
the connection to the Liouvillc model, to the Liouvillc field theory and to Liouville quantum gravity 
measures, is an important direction for further research. In particular, the distribution of the length 
of a segment in Liouville quantum gravity seems to directly connect to our results. 
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Note added: 

After submission, we learned of a recent independent study by D. Ostrovsky [27] who obtained 
a high temperature expansion of arbitrary moments for the [0, 1] problem with no edge charges, and 



conjectured a formula for these moments. Exploiting the integral representation (44 ) for the generalized 
Barnes function Gp(z) together with the following doubling formula: 

G p {2z) = C 2^-^H>n-*G p (z)G fj (z + ^)G p (z + ^)G p (z + ± + | ) (63) 

where Cp is determined from e.g. z = 1, we were able to show that his conjecture is equivalent to our 



formula (47 1. Note however that no discussion of the critical case, duality and freezing is given in [27] . 

We have also shown that using Dirichlet boundary conditions at large distance \x\ = L for the 
2D GFF gives a proper meaning to the problematic Gaussian weight case (it yields a shift 2lnL in the 
second cumulant y 2 , while maintaining all higher cumulants as given in the text). 

Appendix A. the special case of [0, l]_i /2.-1/2 
Here we study the model defined by the partition sum: 

aX a PV(x) 



eP*™ (A.l) 



V(x)V(x') = -2 In \x-x'\ (A.2) 

which is a special case of the interval [0, l] a b defined in the text for a = b = —1/2. We show that it 
corresponds to a REM with a correlation matrix which can be diagonalized in the Fourier basis. 

Using the change of variable x = cos 9, hence dx = sin 9d9 = \Jl — x 2 d9, we see that Z can as well 
be written as an integral over a half-circle, involving a new gaussian random potential with a modified 
correlator: 

Z= dOe^O) (A.3) 



11 



V(0)V(6') = -21n|cos(0) -cos(0')| = ^ - cos(n0) cos(n0') +2 In 2 (A.4) 



n 

n=l 
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where we have used the formula: 



— cos(nA) cos(nB) = — ln(2| cos A — cos B\) 



(A.5) 



n=l 



To define the corresponding REM we now take a grid 9 — 2iri/M for the full circle and take for 
correlation matrix: 



x , 4 

dj = } cos(2wri/M) cos(2n?r j/M) + 2 In 2 



(A.6) 



for i ^ j ■ The sum is still infinite, but it has the nice property that it can be made finite. Indeed using 
that: 



cos(2(fc + mM)m/M) cos(2(£; + mM)nj/M) = cos(2kiri/M) cos(2knj/M) 
we can now rewrite: 



M 



dj = ^2 A » cos(2mri/M) cos(2rmj/M) + 2 In 2 



a„ = - + E( 



1 



1 



m— 1 



n + mM mM 



1 

n 



7 + ^(l + g) 
M 



(A.7) 



(A., 



(A.9) 



One can check numerically that all the eigenvalues are positive. Note that we have subtracted an 
infinite part on the diagonal so now the diagonal element is also well defined: 

M 

Cu = ^A„cos(2mri/Af) 2 + 21n2 (A. 10) 

n=l 

Hence for this particular interval model a = b = — 1/2 we can use the Fourier basis to generate the 
variables on the full circle and take the minimum only for the half circle (i.e. M/2 x M/2 submatrix). 
It remains to be understood how this links to the simplification observed in formula (37 1 in the text, 
and whether there are other examples of such cases where a Fourier basis can be used. 

Appendix B. Some properties of the generalized Barnes function 



Let us first check that the function Gp(x) defined by (44) does indeed satisfy the property (46 1. We 
start from the formula (see 22J 8.341.3, p.! 



00 dt r e~P xt 



t 



1-e 



— + e- t (Px-l)] 



Now, by straightforward algebra (44 1 implies: 
In Gp{x + (3) - In Gp{x) = § In(27r) + 



dt 



*//3) 



-(20x-l) 



Changing now t — > (it, and subtracting (B.l ) gives: 



(x) = lnG fj (x + (3) -]nG p (x) - InT(xP) - ^ ln(27r) 



dt , 



-t/3 



t L l 



+ 0x (e~ tp - e-*) + 



Now using the identity 
dt 

T 



In/3 



(B.l) 
(B.2) 

(B.3) 
(B.4) 

(B.5) 



we see that 



(x) = — /31n/3 which implies 



Statistical Mechanics of Logarithmic REM: Duality, Freezing and Extreme Value Statistics ofl/f Noises generated by Gaust 
with 



Mo) 



dt 



-t0 



t v l 



(B.6) 



In turn, it is easy to see this integral converges and jfi&pi®) = hence </>/3(0) = |ln/3 + </>/3 = i(0). 
Combining altogether we see that (j)f3(x) = (i — (3x) ln/3 + c, where 



dt . e- f 



e 



lim 



die 



, 1 

^e*-l 2 i' 



which using [22] 8.341.1, p.888 yields 



c = lim 

2^0 



1 



lnT(z) +z-(z- 1/2) lnz - - ln27r- 



(B.7) 



(B. 



o * 2 

The last integral is equal to \ In (z/(z + 1 )), and after straightforwardly taking the limit we find finally 
c = — \ ln27r, in full agreement with (46). 



Next we want to obtain the asymptotics. For this it is useful to note that: 



hp{x):=dl\nG p {x)= / ~(e"* - 



dt 



t 



(l-eH*)(l-e-*/0) 



(B.9) 



Exploiting again the identity (B.5l we can rewrite the above formula in a form more convenient 



for applications, see Eq. (52 1 in the text. For example, by changing variables t — t/x in (52 1 we 



immediately can find the asymptotic behaviour for x — > oo at fixed (3 to be given by 

M z)=lnz-^(/?+i 

as long as x 3> max(/3, The same asymptotic behaviour holds for hp(z) in the complex plane for 

5Rz > and \z\ — > oo. 



Note also the useful doubling formula (63 1 for G^x) which we were not able to trace in the 
available literature. 

Appendix C. High temperature expansions 

Appendix C.l. high temperature expansion of REM models 

It is useful to derive high temperature expansions for a Gaussian REM Zm{P) = Y^iLi e~^ Vi with an 



arbitrary correlation matrix ViVj = Cy. One expands: 

Z M ((3) = M - + \f? £ V? - g/? 3 £ V? + 0(/3 4 ) 

iii 

which leads to: 

In Z M { P ) = InM - A £ Vi + ^ _L £ - ^ £ ^ 

- >i E - ^ £ ^ + 2 4 E w») + o(^) 



(C.l) 



(C.2) 



This leads to the average free energy: 



AP 
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and the variance: 

T = ^in 2 z M (p) = ~ E c « + 2 (^2 E(^ + \ C D 



^ E(3C y C ifc + C u C jk ) + ~ E CijCu) + 0(/3 4 ) 

ijk ijkl 



(C.4) 



This result for f 2 is useful to test universality in the sense (ii), i.e. at any temperature. Let us 
examine several cases. 

Consider firs t th e periodic case discussed in the text, where CV, is a cyclic (i.e. circulant) matrix, 
i.e. of the form (57l. Then @ij — Ao/M. Fixing Ao = as we did here, we find that the 

expression for the second cumulant of the free energy simplifies and that it vanishes at (3 = as: 



r- 



2M 2 



TrC 2 + 0(/3 4 ) 



2M 2 



(C.5) 



fe#0 



It is now easy to check that both the discrete circular model (56 1, the sharp model (SM) and the long 
range model (LRM) behave in the limit M — > +oo as: 



(C.6) 



k=l 



i.e. as the continuum circular model for which one has (61 1. This is consistent with the conjecture that 
these models belong to the same universality class at any temperature. Furthermore the coefficient of 
1 can also be obtained, e.g. for the discrete circular ensemble (56), as: 

i i n i i ti i~i . n ri i — - ri . < 

(C.7) 



i teo2Sf5L G «=2y o 2^i ^[21n2|sm(- 



7T" 



Note that the diagonal does not contribute to this limit (its contribution is 0(ln M/M) and that will 
be a general fact. 

For the discrete interval model ( 60 1 we compute the two first terms in the high temperature 
expansion. We see that the terms involving Ca cancel out, as it should. In the limit M — > oo we 
replace remaining sums by integrals and get 



P 



-2 Jo + 2/3 2 



h + 51$ - 61 



0(/3 4 ) , 



where we have defined the integrals: 



dxi / dx2 In \x2 — xi\ 



dxi 



di2 In \x2 — xi | 



I 



dxi 



dx-; 



dx3 ln\x2 — x\ I In \x$ — X\ \ 



(C.8) 

(C.9) 
(CIO) 



Calculation of these integrals give: 

Jo = -3/2 , h 
which gives the final result: 

F = 3 + /? 2 



7 / _17_tt 2 _ 
2' 6 18 



2 2 9 

-7T 



0(/3 4 ) , 



(C.ll) 



As we show below this coincides with our analytical prediction from the continuum model, see ( C 17 1 
below. 

As we see in Fig. [IT] for the discrete interval model (60 1 at finite M f 2 (f3 = 0) is smaller than 
3. In fact, one can add a Wm on the diagonal in (60 1 so as to tune this value to exactly 3 for any 
M, without changing the universality class (i.e. Wm goes to zero fast enough). One could try to 
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systematize this idea, e.g. to add to the correlation matrix of the discrete model some other matrix, 
subdominant in the limit M — ► oo, so as to fit the lowest orders coefficients in (3 V to their actual value 



for the continuum model - those are given below for the interval, sec formula (C.17 1. We have checked 
for the circular case that it can be easily implemented up to p — 2. Whether this will allow to select 
better discrete models with faster convergence even at lower temperature is left for future studies. 

Concerning the class of model ([5]), we can similarly check that for the associated discrete REM, 
i.e. Cij = -21n|2sin(§/(0i) - \f{0j))\ for i ^ j and 9. t = 2ni/M one has: 

f27T ja p2-rr 

lim -^ V C„ = -2 



dO l 
2tT 



— ln[2|sm( 



1 , 
-a' 
2 



+ 0(a 4 ) (C.12) 



for f(6) = 8 + asin(#), hence at small a, p C = |a 2 + 0(a 4 ) + 0(/3 2 ) as announced in the text, and 
there is no universality valid at all temperature (the universality class in the sense (ii) defined above 
depends on the function f(6)). 

Finally the same expansion (C.4) holds for any continuum REM of the form Z = J dxp{x)e~^ vtyX ^ 
and can be obtained from the above just replacing i n ~* (J dxp(x)) n Xci x an< ^ ^*i>»2 by its 

continuum expression Cixx^x^)- 



Appendix C.2. high temperature expansion of the analytical result for the interval 



Let us derive the high temperature expansion of our analytical result (50 1- (52 1 for the second cumulant 



of the free energy. Since an independent method also exists to obtain this expansion, as displayed in 
the discussion above, this constitutes a check of our solution in the high temperature phase. For this 
we need to use (52 1 for x — ai[3+ct2jj, where ai.2 are given positive constants. By introducing r = t//3 
we have: 



ai fl- 



an 



= In 



a 2 



ln[l + ^/? 2 
a 2 



dr 



-a 2 T-a 1 /3 t 



1 - 



(l_ e -/3 2 r)( 1 _ e -r) 



(C.13) 



which can be easily used to expand in powers of (3 2 . In particular, the leading term from (50 1 is a 
constant given by: 

dr 



r- 



21n3 + ln2- 



1 



1 



2e-*'' 



where the integral can be computed by pieces using the formulas (see 3.311.7 [2"2~]): 



dr 



1 



ip{v) — V ; (m) 5 an d 



dr 



(C.14) 



In- (C.15) 



combining altogether we find: 

P\(3 = 0) =41n2 - [^(1) + V(2) - 2^(3/2)] 



(C.16) 



in agreement with the result obtained above in (C.9) by a direct method 



This expansion can be carried to higher order. Using mathematica and some heuristics we find 
that it can be put in the form: 



P C = 3 + ( ^ 2 - I )P 2 + E(- 1 ) fc+l3fc (C(* +!)-(!+ T 

'* k=2 " 



2/,- 



(C.17) 



where the are the Bernouilli numbers (B^ — for k odd). 
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